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Q , Abstract 

\o . 

■ We investigate the leading power corrections to the decay rates and distributions in 

, the decay B Xsi~^£~ in the standard model (SM) using heavy quark expansion (HQE) 

CL|| in (l/nib) and a phenomenological model implementing the Fermi motion effects of the 

Q^' b quark bound in the B hadron. In the HQE method, we find that including the lead- 

D . ing power corrections the decay width T(B — > Xsi^i~) decreases by about 4% and the 

branching ratio B{B Xsi~^i~) by about 1.5% from their (respective) parton model 
values. The dilepton invariant mass spectrum is found to be stable against power cor- 
^ ■ rections over a good part of this spectrum. However, near the high-mass end-point this 

distribution becomes negative with the current value of the non-perturbative parameter 
A2 (the Ai-dependent corrections are found to be inoccouos), implying the breakdown of 
the HQE method in this region. Our results are at variance with the existing ones in the 
literature in both the decay rate and the invariant dilepton mass distribution calculated 
in the HQE approach. As an alternative, we implement the non-perturbative effects in 
the decay B Xsi^i~ using a phenomenologically motivated Gaussian Fermi motion 
model. We find small corrections to the branching ratio, but the non-perturbative ef- 
fects are perceptible in both the dilepton mass distribution and the Forward-Backward 
asymmetry in the high dilepton mass region. Using this model for estimating the non- 
perturbative effects, modeling the dominant long distance (LD) contributions from the 
decays B ^ Xg + {J/ipjtp' , ...) Xsi~^i~, and taking into account the next-to-leading 
order perturbative QCD corrections in 6 ^ s£~^i~ , we present the decay rates and distri- 
butions for the inclusive process B Xsi~^i~ in the SM. 
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1 Introduction 



Rare B decays B Xsi~^i~ and B — > Xg'y are well suited to test the SM and search for physics 
beyond the SM. In the SM, such processes are governed by the Glashow-Iliopoulos-Maiani 
(GIM) mechanism and their rates and distributions are sensitive to the top quark mass and 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements P|. First measurements of the decay 
rates for the exclusive decay B K* + 7 @ and the inclusive decay B Xg'j have been 
reported by the CLEO collaboration. At the partonic level, the complete leading order (LO) 
anomalous dimension matrix involving the b —>■ s'j decay was calculated in in the context 

of an effective five-quark theory. First calculations of the gluon bremsstrahlung and virtual 
corrections, which are part of the next-to-leading-order perturbative QCD improvements, were 
reported in [§]-[|10| (see also |Tl], 0). The NLO virtual corrections to the matrix elements have 



been completed in |13|. A first calculation of the hitherto missing NLO anomalous dimension 



matrix has been recently reported in |T^. Leading power corrections in (l/mj,) to the decay 
rate r{B ^ Xg'y ) (and r{B ^ Xiui), which is often used to estimate the branching ratio 
B{B Xs'~f )) have also been calculated in the heavy quark expansion (HQE) approach [|15|, p^. 



A quantitative measure of the rapport between experiment and present estimates in the SM is 
the CKM matrix element ratio |V(s|/|Vc6| for which a value |Vts|/|V"c6| = 0.85 ± 0.12 (expt) ± 
0.10 (th) has been obtained from the inclusive decay rate for B Xg^ [|l^, in agreement with 



the bounds obtained from unitarity of the CKM matrix |18 |. 

It is known that the inclusive energy spectra in the decays B Xiu^ and B —>■ Xgj 
are not entirely calculable in the HQE framework [ll6|,ll^-[^. In particular, the end-point 



energy spectra are problematic in that the energy released for the light quark system in the 
decay b qX (here X = 7 or a dilepton pair) is not of order rrih but of order A, where A = 
mB—TTLh = O(Aqqj)). Hence, the expansion parameter in the HQE approach, which is formally 
of 0(1/(5^) = 0{l/ml), near the end-point gets replaced by a quantity which is of 0{1/K'^) = 
0{1/Amf,) ^ 0{l/ml), implying the onset of the breakdown of the HQE method. To make 
contact with experiments one has to smear the energy spectrum in question over an energy 
interval sufficiently larger than Aqqj). Thus, direct comparison of theoretical distributions 
with experiments requires additional input in terms of phenomenological models, e.g., the 



Gaussian Fermi- motion model of [^, which incorporate such smearing. The smearing effects 



are very important in B ^ Xg'-f ^ |T2[ and not negligible in the lepton and hadron-energy 
spectra in the decays B Xc{Xy)ii^e, either |^-||2^. Alternatively, one may have to resort 
to a resummation of the power corrections near the end-point [ll9| , p5| . Such resummations, 
however, remain so far inconclusive. 

In this paper, we address the related FCNC process B Xsi~^i~, i = e,fi. (Since we 
neglect the lepton masses in our calculation, our results are not applicable to the decay B — >■ 
XsT~^T~.) The SM-based short distance (SD) contribution to the decay rate for the partonic 
process b si'^i~, calculated in the free quark decay approximation, has been known in the 



LO approximation for some time p6|. In the meanwhile, also the NLO perturbative QCD 
corrections have been calculated which reduce the scheme-dependence of the LO effects in 
these decays [^, In addition, long distance (LD) effects, which are expected to be very 



important in the decay B Xgi'^i [2^, have been estimated from data on the premise that 



they arise dominantly from the charmonium resonances J/ip and ip' through the decay chains 
B — > XsJ/iplip') — > Xgi^i^ . Higher resonances {tp", ip'", ...) also contribute though at a reduced 
level. Estimates of the LD effects away from the resonance regions involve specific assumptions 
about the g^-dependence of the relevant vertices, which at present can only be obtained in 
specific models - p2|. 



The particular aspect we are interested in is an estimate of the non-perturbative effects on 
the decay distributions in B ^ Xgi^i' , which take into account the i?-hadron wave function 
effects and incorporate the physical threshold in the final state on the underlying partonic 
calculations. This effects both the SD- and LD-contributions, and to the best of our knowledge 
has not yet been calculated. Closely related to this aspect is the question of power corrections 
to the parton model decay rates and spectra which have been calculated for the SD-part of the 
dilepton invariant mass distribution in B ^ Xsi'^i" by Falk, Luke and Savage (henceforth 
this paper is referred to as FLS) using the HQE approach. We reevaluate these corrections in 
this paper, reaching different results and conclusions than in the FLS paper which we specify 
later. 

From the power corrections calculated in the HQE approach for the decays B Xg'y and 



B Xiui |T^, [T^, we recall that there are no leading, i.e., 0(l/mfe), corrections in the inclusive 
rates. Likewise, in the decay B Xsi'^i~ , the first non- vanishing corrections to the inclusive 
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rates are of 0(l/m^). Furthermore, the dilepton mass spectrum in B ^ Xgl^l is found to be 



well behaved in the HQE framework in the entire dilepton mass range in FLS . In particular, 
the high dilepton invariant mass spectrum in the parton model is found to receive moderate 
power corrections, typically 0(10%), increasing the dilepton yield in S — > Xs^^^~ (see Fig. 2 
in pH]). This result differs qualitatively from analogous power corrections in the lepton energy 



spectra in B ^ Xiu^, which are large and negative near the end-points (see, for example. 



Figs. 5 - 8 in the paper by Manohar and Wise ||T6|). In addition, taking the V — A limit in the 
matrix element for B X^i^i', the differential distributions and decay rate in this process 
can be related to the corresponding quantities in the semileptonic decay B Xlv^. The power 
corrections in the latter decays have been calculated and discussed at great length by Bigi et 
al. in |T5[ and by Manohar and Wise [Q. We are of the opinion that both the power corrected 
dilepton spectrum and the inclusive decay rate V{B ^ Xs£^i~) obtained by integrating this 



spectrum in FLS are at variance with the results in [^, [T^ in this limit (see Appendix C). 
In view of the impending interest in the decay B Xgi'^i", in particular the dilepton mass 
spectrum and the Forward-backward asymmetry involving £^ versus i~ [Q, which have been 
put forward as precision test of the SM in the FCNC sector and hence a possible place for 
discovering new physics |35|, |3^ , we have recalculated the power corrections in this process in 
the SM using the HQE approach. 

To that end, we have computed the Dalitz distribution, d^B/dsdu, for the decay B — >• 
Xsi~^i~ (see section 2 for the definition of these variables), taking into account the NLO per- 
turbative QCD correction in and the leading l/m^ corrections in the HQE approach. In 
doing this, we have also kept the s-quark mass effects. Integrating over one of the variables, 
the resulting expressions for the dilepton invariant mass and the FB asymmetry are derived. 
While the power-corrected FB asymmetry in S — > Xgi^i^ is a new result, not presented earlier, 
our expression for the power-corrected dilepton mass distribution is not in agreement with the 
one presented in FLS [Q. Since the derivations of the final results for dr{B Xgi'^i" )/ds 
and the FB asymmetry A{s) are rather involved, we have decided to give the details of the 
calculations so that they can be checked stepwise and the source of this discrepancy pinned 
down accordingly. Some checks of our results in the limiting case mentioned above are already 
possible and have been carried out. In particular, we are able to derive the results in O, IIH] 
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taking the appropriate limit of our expressions in 5 — > Xs^"''^~(see again Appendix c). 

We find that the final-state distributions in i? ^ X^f^f" are not calculable entirely in the 
HQE approach, as the dilepton mass distribution becomes negative in the end-point region. 
While this defect may be resuscitated by resummation of the HQE-power corrections, we do 
not attempt this here. Instead, we estimate the non-perturbative effects on the decay rates and 
distributions in B —>■ Xsi^i~ by invoking the Gaussian Fermi motion model [^. This model 
has been used successfully in the analysis of the lepton energy spectrum in the semileptonic 
decays B Xiug and the photon energy spectrum in B ^ X^j As pointed out in 
]16| on the example of B Xiui, this model reproduces the effect due to the kinetic energy 



term Ai in the HQE approach, if the 6-quark mass is appropriately defined, but there is no 
analogue of A2 (the matrix element of the magnetic moment operator) in the Fermi motion 
model. The distributions in the two approaches (HQE and the Fermi motion model) are hence, 
in general, different, which is most noticeable near the end-points. By construction, there are 
no negative probabilities encountered in the Fermi motion model and the final state thresholds 
can be correctly incorporated. 

This paper is organized as follows: In section 2, we derive the double differential distribution 
d'^B/dsdu for the decay B X^f+f", including the explicit 0{as) and the leading power 
corrections in 1/ rrib, giving in Appendix A the individual contributions to the structure functions 
from several contributing sources governing these decays. Some of the lengthy expressions 
obtained in the derivation of the HQE-improved Dalitz distribution are displayed in Appendix 
B. The power-corrected dilepton invariant mass distribution and the FB asymmetry in B 
Xgi^i^, together with their simplified versions in the limit = 0, are also given in this section. 
We also present here numerical comparisons in the two quantities of interest between the parton 
model and the HQE-approaches, as well as differences between our result and the one in FLS 



33| . In Appendix C, we present the limiting case of our results for B Xsi^i and compare 



them with the existing ones in the literature [|T5|, [T^. In Appendix D, we show (a peripheral 
result) that the energy asymmetry defined in |]5B| and the FB asymmetry introduced in 
are related. In section 3, we implement the B- meson wave function effects and the physical 
threshold on the final state in B Xs£^£~, using the NLO-corrected parton distributions and 



the Gaussian Fermi motion model |22|. Since the calculation of the FB asymmetry in this model 
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involves some non-trivial kinematic transformations, we have given the details in Appendix E. 
The LD-contributions in B ^ Xgi'^i" are estimated in section 4, using data on vector meson 
intermediate states B ^ V + Xs, where V = {J/ip{lS), ...,ip{6S)). The resulting dilepton 
mass spectrum and the FB asymmetry, including the wave-function and LD-effects, are also 
presented here. We conclude with a discussion of our results and possible improvements of the 
LD-effects in 5 — > Xgi'^i' in section 5. 

2 Power corrections to the dilepton invariant mass dis- 
tribution and FB asymmetry in B ^ Xs£^£~ 

We start by defining the various kinematic variables in the decay b{pb) —>■ s{ps)+i'^ ip-\-)+i~ (p-) ■ 

u = -{Pb - P+f + {Ph - P-f, 

u{s,ms) = \J{s - {rrib + msY) {s - {rrib - msY). (1) 

For subsequent use, wc note that p± = {E±,p±), and = is the momentum transfer 

to the lepton-anti-lepton pair (hence = s). We also define the four- velocity of the h quark, 
■^/i — iPb)fM/nT'b, which we shall take subsequently to be the same as that of the B hadron, 
i^n — {pb)i^/Mb- Finally, we introduce the scaled variables s and u 

s 

S — 2 ) 

mi 

XL 

^ = -2 =2^-(P+-P-)' (2) 
which in the decay b —>■ si~^i~ are bounded as follows, 

—u{s, rhs) < u < +u{s, rhs) , 
u(s,ms) = ^J[s - (1 + ms)2] [s - (1 - ms)^] , 

Am^ < s < {l-rhsf . (3) 

where rhi andp,^ are the scaled quark masses and scaled momenta, respectively, rhi = mi/mh, Pi — 
Pi/rrib. 
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2.1 NLO-corrected amplitude for b st^l in the effective Hamil- 
tonian approach 

Next, the explicit expressions for the matrix element and (partial) branching ratios in the 
decays h si^i~ are presented in terms of the Wilson coefficients of the effective Hamiltonian 
obtained by integrating out the top quark and the bosons, 

AC 

Heffib -^s + X)=neff{b^s + ^)~ ^ytlVtb [C9O9 + CloOlo] , (4) 

where 

n^ffib . + 7) = -^v*v,,j:c^ii^mf^) ■ (5) 

Here, Vij are the CKM matrix elements and the CKM unitarity has been used in factoring out 
the product V^lVtb- The operator basis is chosen to be (here fj, and z/ are Lorentz indices and a 
and P are colour indices) 

01 = {sLalfMbLa){cL(il''CLp), (6) 

02 = {sLa'yfibLI3){cLl3l^CLa), (7) 

03 = {sLa^t^bLa) i^Lp^^qifd) , (8) 

q=u,d,s,c,b 

Oi = {sLali^bLp) Yl iQLfsYqLa), (9) 
q=u,d,s,c,b 

O5 = (sLalt^bla) J2 (QRIbYQRp), (10) 
q=u,d,s,c,b 

Oq = {sLalt^bLo) J2 (^i?/37^9i?a), (H) 

q=u,d,s,c,b 

O7 = Y|^s,fT^,(mbi? + m,L)6„F^^ (12) 
Os = j^s^T^pa^,{m,R + msL)bf,G''^'\ (13) 

(14) 

where L and R denote chiral projections, L{R) = 1/2(1 =F75)5 and the two additional operators 
involving the dileptons are: 

= T^Sal''LbJ-f/, 
lOTT^ 

^10 = -i^^cYLbJj^^.i . (15) 



The Wilson coefficients are given in the literature (see, for example, |28[|). 
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With the help of the effective Hamiltonian in eq. the matrix element for the decay 
b —>■ si~^i~ can be written as, 



M{b^si+i- ) 



V27r 



-2Cf (sicr^, ^i^sL + rrihR) bj [h^^ I 



(16) 



We have kept the s-quark mass term in the matrix element explicitly and this will be kept 
consistently in the calculation of power corrections and phase space. The above matrix element 
can be written in a compact form, 



M{b^ sfM- ) = ^ V*V,, (t\ L^" + T\ L^') 



(17) 



with 



= h,L{R)L 



b. 



(18) 
(19) 



We recall that the coefficient Cg in LO is scheme- dependent. However, this is compensated by 
an additional scheme-dependent part in the (one loop) matrix element of O9 [p7| , p8[ . We call 
the sum Cg , which is scheme-independent and enters in the physical decay amplitude given 
above, with 



cfis)^CMs)+Y{s). 
The function Y{s) is the one-loop matrix element of Og and is defined as 



(20) 



~gil,s) (4^3 + 4^4 + 3^5 + ^6) 



vis) 



g{0,s) (C3 + 3C4) 



+ -(3C3 + C4 + 3C5 + C6) 

y 

-4(3Ci + C2-C3-3C4) 
1 + ^ujis) . 

TT 



(21) 
(22) 



Here, ^ is dependent on the dimensional regularization scheme with 



(NDR), 
-1 (HV), 



(23) 



in the naive dimensional regularization (NDR) and the 't Hooft-Veltman (HV) schemes. The 
function uj{s) represents the 0{as) correction from the one-gluon exchange in the matrix element 

of Og 



2 4 
2s(l + s)(l -2s) 



- In s ln(l — s) — — ln(l 

3 ^ ' 3(1 + 2s) ^ 



Ins + 



5 + 9s - 6s2 



3(1 - s)2(l + 2s) 6(1 - s)(l + 2s) ■ 

The function g{z,s) includes the charm quark-antiquark pair contribution 



(24) 



28 



X 



e(l-y)(ln 



9(0, §) 



l + y/T^ 

_8_ _ 8 nib, 
27 9 ^ V ' 



in) + Q{y — 1)2 arctan 



4, . 4. 
- m s H — tir 
9 9 



(25) 



(26) 



where y = Az^ /s. With the help of the above expressions, the differential decay width becomes. 



dr 



1 G/a' ,2 ci3p+ dV 



\v:sVtb\ 



{w 



fJ,U 



(27) 



2Mb 27r2 ' (27r)32^+ (27r)32E_ 

where W^i, and L^^, are the hadronic and leptonic tensors, respectively. The hadronic tensor 
W^/^ is related to the discontinuity in the forward scattering amplitude, denoted by T^J^, 
through the relation W^i, = 2ImT^i,. Transforming the integration variables to s, u and v ■ q, 
one can express the Dalitz distribution in 6 — > si~^i~ (neglecting the lepton masses) as: 
dr 1 Gp'^a'^ mb^ 



with 



dwds 2mb 2 7r2 256 tt^ 



rpL/R 



\V:M' 2 Im / d(.; . q) [t\^ L^'" + L^'" 



d'ye~^^-y (b |t {rij/^(y), r2^/^(o)}| b) , 

E [v''^'^{p+)Yu''^'^{p-)] [M^/^(p-)7"t;^/^(p+) 

spin 

2 - g^^^'{p+■p.)T^e^^""'p+^p 



(28) 



(29) 



(30) 



where T^f^^^^ = TaJ/-^ = T^/^, given in eq. The Dalitz distribution eq. (|8|) contains the 

explicit 0(as) -improvement, and the two distributions in which we are principally interested 
in can be obtained by straight-forward integrations. 
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2.2 Leading power (l/m^,) corrections in the decay B X^t^l 

The next task is to expand the forward scattering amphtude T^j^ in the inverse powers in l/mf,. 
Suppressing the Lorentz indices for the time being, this expansion can be formally represented 
as: 



1 

rUb 



{B\Oo\B) + -^{B\Oi\B) 
2 rrib 



{B\02\B) + -- 



Arab" 



(31) 



and the expressions for the operators Oq^Oi and O2 are given explicitly in ^3[. They are 
obtained by expanding the propagators in the Feynman diagrams contributing to the time- 
ordered product on the l.h.s. of the above equation (see Fig. 1 in |^), using pb^l = ''nfihV^ + fc^, 
fixing the four-velocity of the external h quark field to be f ^ and treating the components of 
the "residual momentum" to be much smaller than m^. 

As is well known, the leading power corrections can be parametrized in terms of the matrix 
elements of the kinetic energy and magnetic moment operators, called Ai and A2, respectively, 
and defined as, 



B 



hiiDYh 



B 



B) = 2MbXi 
B) = 6Mb\2 



(32) 



where B denotes the pseudoscalar B meson, is the covariant derivative and G^j^ is the QCD 
field strength tensor. The two-component effective field in the HQE approach h{y) is related 
to the QCD field b{y) through the expansion. 



2mb 



h{y), 



(33) 



where /D = D^'-y^. The parameters Ai and A2 are related through the quantity A to the hadron 
masses [^ , 

-r Ai + 3A2 

rriB = rrib + A h 

Irrib 

- Ai — A2 

rriB* = rrib + A h .... 

2mh 



(34) 



^From the B — B* mass difference, one obtains A2 — 0.12 GeV^. The quantity Ai has been 



determined from QCD sum rules ^ and data Its present value is subject to a certain 
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theoretical dispersion, estimated somewhere between Ai = — (0.52±0.12) GeV^ (Ball and Braun 
in to Ai = -(0.10 ± 0.05) GeV^ (Neubert ||). 

Concerning the definitions of the operators in eq. (|3T|) , we follow the prescription given in 
p3| , in which the leading operator Oq is defined in terms of the "full" four- component field 

Ky), 

Oo{y) = hr,{/j- j^ + ms)T2b, (35) 

where x = l + s~2{v-q)— ml + ie. The other two subleading operators Oi and O2 are, 
however, written in terms of the two-component effective fields h{y), which is related to the 
field b{y) through the expansion given in eq. (|33|). Of these, the expression for O2 involving 
the expansion of the one-gluon graph is obtained by a non-trivial derivation, which we have 
checked, and it agrees with the one given in eq. (3.8) of [Q. (Likewise, we agree with the 
expression for Oi given in eq. (3.6) of [@.) For the sake of completeness, we give below the 
explicit expression for Oi and O2 given in 



Oi{y) = hva^T^tD^h -^{v~ qrhT,{/j- ^ + m,)V2%DJi, (36) 
X x^ 

and 

16 - 4 - 

02{y) = -^{v - qT{v - qfhT^{^- ^ + m,)V2iD^iDph - -^hT^{^- ^ + rn,)T2{iDfh 

- \{v-qfhTi-i''T2{iDJDp + iDpiDo)h + \mshTiio^pT2G'''^h (37) 

3y Ob 

+ ^le^^'^^iv - q)xhTr^,^^T2G^^h + 'h{^^Ta''T2 + T ^^^T 2l'')iD piD 

- - g)"/^7^ri(/;- ^ + m,)T2iDpiD^h - ^{v - qThT^{^- ^ + m,)T2l''iD^iDph. 

Using Lorentz decomposition, the tensor T^y can be expanded in terms of three structure 
functions, 

T^u = -Ti g^u + T2 Vf, Vu + T3 ie^uap v°' f , (38) 

where the structure functions which do not contribute to the amplitude in the limit of massless 
leptons have been neglected. After contracting the hadronic and leptonic tensors, one finds 

T^/ V L'^^''"^ =m^l^2s Ti^/^ + [{v ■ qf - ^u' - T//^ t § u Tg^/^} . (39) 

We remark here that the T3 term will contribute to the FB asymmetry but not to the branching 
ratio or the dilepton invariant mass spectrum in the decay B ^ Xg 
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The results of the power corrections to the structure functions Tj can be decomposed into 
the sum of various terms, denoted by t/-'^ which can be traced back to well defined pieces in 
the evaluation of the time-ordered product given above: 

T,iv.q,s)= (40) 

j=0,l,2,s,g,5 

The expressions for t/^'^ {v.q,s) calculated up to 0{MB/ml) are given in Appendix A. They 
contain the parton model expressions Ti^\v.q, s) and the power correction in the HQE approach 
which depend on the two HQET-specific parameters Ai and A2 defined in eqs. (|32|) . Note that 
the s-quark mass terms are explicitly kept in T^^\v.q, s). 

From the expressions for t/^'^ given in Appendix A, we see that Tf\i = 1,2,3) are of order 
MB/{mb) and the rest t/^^ , t/''\ t/^\ t/'^ and t/^^ are all of order MsAi/mb^ or MsAa/mfc^. 
Since the ratio M^/mfe = 1 + 0(l/m;,), we note that the Dalitz distribution in i? — > XgC'^Cr has 
linear corrections in l/rrib- The origin of the various terms in the expansion given in eq. ( ^OD is 
as follows: 

• The contributions to t/^"* come from the matrix element of those terms in the operator 
O2, which originate from expanding the spinor of the heavy quark field b{x) in terms of 
the spinor of the heavy quark effective theory h{x). 



The remaining contributions from the matrix element of the operator O2 are denoted by 

ina1 

(2) 

diagram and the rest being T{ . 



and t/^\ with t}^^ originating from the matrix element of the one gluon emission 



• The contributions denoted by t/^'' arise from the matrix element of the operator Oi. In 
the leading order in {l/rrib) this matrix element vanishes, but in the sub-leading order 
it receives a non-trivial contribution which can be calculated by using the equation of 
motion. 

• Finally, the contributions t/'*'' arise from the matrix element of the scalar operator bb. 
Concerning the last point noted above, we recall that the scalar current can be written in terms 



of the vector current plus higher dimensional operators as |T^ 



bb = v^b^b + -^h UiDY - {v.iDf + s^^'G^, h + (41) 
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with G^^ = [iD^, iD^] and s^^ = {-i/2)a^"'. We note that in deriving t/'\ use has been made 
of the conservation of the 6-number current in QCD, which yields the normahzation: 



{B\h,b\B) = 2{pb] 



(42) 



Finally, after doing the integration on the complex plane v ■ q (see Fig. 1 in for the 
analytic structure of T^jy and the contour of integration), we derive the double differential 
branching ratio in B ^ Xgi'^i" ■ The result can be expressed as, 
dB 



ds du 



1 - -2 -2 

1 — m„ — s — u 



(2 Ai(-1 + 2m^ - - 2 s + s^] 



+3 A2(-l + 6m^ - 5mt - 8 s + 5 s 



+ 



4(1 



^ 4 I ^ 



N^eff|2 I 1 2 



-2-2 I -2 I -2-2 



-- (2 Xi{-1 + ml + m^^ - ml + 2 s + lOm^s + + m^s^ 



2s2^ 



+3 A2(3 + 5mf - 3m'; - 5m° + 4 s + 28mf s + 5 + 5?fif s 



efri2 



s(l + m^) - (1 - mlYj + 2Ai(-l + 2m^ - + s + m^s) 

.effN r^eff 



+A2(5mf - 5m^ + 2s + 5mfs) Re{C^'') 



+2 



2 + \i + 5\2\usRe{cf) Cio 



+4 [2 (1 + ml) + Ai(l + ml) + A2(3 + 5m^) 



u 



ReiCw) Cf^ 9 



—Ei{s,u)6 u{s,ms)'^ — u^ 



Eois, ii,) 5' 



u[s,ms 



u 



- \2 -2" 

u{s^ms) ~u 

(43) 



where Ai = Ai/m^ and A2 = \2lTnl. The auxiliary functions Ei{s,u) {i = 1,2), introduced here 
for ease of writing, are given explicitly in Appendix B. The boundary of the Dalitz distribution 
is as usual determined by the argument of the ^-function and in the {u, s)-plane it has been 



specified earlier. The analytic form of the result (^31) is very similar to the corresponding 
double differential distributions derived by Manohar and Wise in [16| for the semileptonic 
decays B —>■ {Xc, Xu)ii^e- Further comparisons with this work in the V — A limit for the single 
differential and integrated rates are given in Appendix C. 

It has become customary to express the branching ratio for B Xsi~^i~ in terms of the 
well-measured semileptonic branching ratio Bsi for the decays B (Xc,X„)£z/^. This fixes the 
normalization constant Bq to be. 



Bn = B 



'si 



3«2 




1 


167r2 




2 


f{mc)n{m^) 



(44) 
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where 



/(mj = l-8m^ + 8m^-m^-24m^ Inrhc (45) 
is the phase space function for r{B — ^ XJu) in the lowest order (i.e., parton model) and the 
function n{rhc) accounts for both the 0(0;^) QCD correction to the semi-leptonic decay width 
22| , p7| and the leading order (l/m?,)^ power correction |jl5|. Written explicitly, it reads as: 

(46) 



2as{mb) , . - h{m 
1 9{mc) + 



2m 



2 ' 



where the two functions are: 



9{mc 
h{rhc 



(-^-|)(l--e)^ + ^ 



Ai + 



A, 



-9 + 2Arhi - 72mt + 72m° - 15m, - 72m^ Inm, 



(47) 



Finally, after integrating over the variable u, we derive the differential branching ratio in 
the scaled dilepton invariant mass for B Xsi~^i~ , 



as 



+m^s + m^s — rhgS — 3s — 2m^s — 3m^s + 5,s +5171^3 —2s 



Ai 



u[s,ms) 



+ ( 1 — 8ml + 18m^ — 16m^ + 5m^ — s — 3m^s + 9m^s — 5m^s — 15,s^ — ISm^s^ 



-15m^s^ + 255^ + 25mls^ - lOs^ 



A, 



cf\^ + \a 
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+ 



^m(s, m,)(2(l + m^)(l - m^)^ - (1 + 14m^ + m^)s - (1 + ?fi^)s^^ 



+ -(2 - 6m^ + 4m^ + 4m^ - 6m« + 2ml^ - 5s - 12mp + 34m^s - 12m^s - 5m*s + 3s' 



Ai 



u(s,m^J 



+29mf s^ + 29m^s^ + 3m°r + s'^ - lOmf s^ + m^s'^ - s* - s 
+20m^ - 12m^ - 14m^ + 10m^° + 3s + 16m^s + 62mfs - 56m^s - 25m«s + 3s 



+ 4 -6 + 2m^ 



+73mp^ + lOlmts' + 15m^s' + 5s^ - 26m^s'' + 5m^s^ - 5s* - 5mls* 



+ 



s,^ 0,2 

A2 

u{s,ms) 



\C'- 



eff 12 



u{s, ms)((l — mlY — (1 + 'nT''i)s) + 4(1 — 2ml + — s — tt^'Is) u{s, rhs) Xi 
+4 (-5 + 30m^ - 40m^ + 15m^ - s + 21m^s + 25m^s - 45m^s + 13s^ + 22m^s^ 

A2 



+45mts^ - 7s^ - 15m^s^ 



u[s, mg 



Reicf)cf 



(48) 



Another interesting quantity is the FB asymmetry defined in [0, |35 

d^(s) /•! d^B 



ds 



ds dz 



dz 



d^B 
-1 dsdz ' 



(49) 
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where z = cos6 is the angle of measured w.r.t. the 6-quark direction in the dilepton cm. 
system. The leading power corrected expression for the FB-asymmetry A{s) is: 

dA{s) 



ds 



-2 Bo 



2{u{s, ms)yrs + -(3 - 6mf + 3m^ + 2s - 6m^s + 3r)Ai 



+s (-9 - Gm'i + 15m^ - 14s - 30m^s + 15s') Aa Re{C^^) C 
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+ 



4(n(s, ms)Y(l + ) + - (1 + m^) (3 - 6mf + 3m^ + 2s - 6m^s + 3r )Ai 



-2(-7 - 3m' - 5m^ + 15m^ - 10s - 24m's - 30m^s + 9s' + 15m's') As 



e 



/^From the experimental point of view, a more useful quantity is the normalized FB-asymmetry, 
obtained by normalizing dA/ds with the dilepton mass distribution, dB/ds, 



dA _ dA dB 

ds ds ds 



(51) 



This asymmetry, which we recall is defined in the dilepton c.m.s. frame, is identical to the energy 



asymmetry introduced in [Q, which is defined in the B rest frame, as shown in Appendix D. 

The results derived for the 0(as)-improved and power-corrected Dalitz distribution, dilepton 
invariant mass, and FB-asymmetry in B ^ Xgi'^i' are the principal new results in this section. 
It is useful to write the corresponding expressions in the limit = 0. For the dilepton invariant 
mass distribution, we get 



-rr = 2Bo 
ds 



-(l-s)'(l + 2s) (2 + Ai) + (l 



15s' + lOs^ ) A. 



r^eS\2 I \ry |2 
*-"9 + *-^10 



+ 
+ 



\ci 



eff 12 



-(1 - s)'(2 + s) (2 + Ai) + 4 (-6 - 3s + 5s=^) A2 
4(1 - s)'(2 + Ai) + 4 (-5 - 6s + 7s') aJ Reicf) cf 



(52) 



The (unnormalized) FB asymmetry reads as. 



2(1 - s)'s + -(3 + 2s + 3s')Ai + s (-9 - 14s + 15s') A2 



Reicf)C,o 



+ 



4(1 - s)' + -(3 + 2s + 3s')Ai + 2(-7 - 10s + 9s') A2 
3 



ReiCo) cf 



(53) 



A direct comparison of our result for the dilepton invariant mass distribution given in 
eq. (I52D above can now be made with the differential decay width dr(i? —>■ Xsi^i~ )/ds derived 



in eq. (3.21) of the paper by FLS ||3^. To that end, one has to take into account the (obvious) 
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normalization difference between the decay width and branching ratio, rewrite the quantities 
and used in FLS |Q in terms of the Wilson coefficients Cf , Cg and Cio used by us, 
with A^^^ = Cg^ ± Cio and B^^^ = —2Cf^, and drop the explicit 0(as)-improvement in the 

rr 

coefficient Cg , as FLS did not include it in their calculations. The resulting expression is: 



dB 



ds 



-(1-5)(1 + 2S) 



+ ^(5 + 3s - 2s2)Ai + ^(1 + 15s - 10s' ) A2 

D 2 



c 



efri2 



IC'ioP 



+ 
+ 



;i + s)Ai - 10(1 + s)A2 



\c 



efri2 



4(1 - s) - 2(-- + s)Ai + 2(5 - 7s)A2 



Re{cf) cf 



We would like to make the following observations: 



(54) 



The results derived here (eq. (52)) and in FLS |33| (eq. (54)) reproduce the known 
parton model expression for the dilepton invariant mass distribution in the limit Ai — > 
and A2 0. 

The power corrections themselves, i.e. the expressions multiplying the constants Ai and 
A2, are different in the two derivations. 

The power-corrected dilepton invariant mass distribution derived by us retains the char- 
acteristic 1/s behaviour following from the one-photon exchange in the parton model, in 
contradiction to the observations made in |Q . This difference can be seen by comparing 
the two expressions multiplying the Wilson coefficient jC^^p. 

Leading order power corrections in the dilepton mass distribution are found to be small 
over a good part of the dilepton mass s. However, we find that the power corrections 
become increasingly large and negative as one approaches s gmax ^ Since the parton 
model spectrum falls steeply near the end-point s §"^-0.^^ f^j^i^ leads to the uncomfortable 
result that the power corrected dilepton mass distribution becomes negative for the high 
dilepton masses - in contradiction to the observations made in ||33|. We show in Fig. ^ 
this distribution in the parton model and the HQE approach, using the central values of 
the parameters in Table ffl. 



15 



We note that the correction proportional to the kinetic energy term Ai renormahzes the 
parton model invariant mass distribution multiplicatively by the factor (1 + Ai/(2m^)), 
i.e. no new functional dependence in s is introduced (moreover, this factor is hardly 
different from 1). Hence, the negative probability near the end-point is largely driven by 
the magnetic moment term A2. 



• A comparison of the dilepton mass spectrum resulting from eq. (|52D of this work and 
eq. (3.21) in FLS |Q (i.e. eq. (54) given above) is shown in Fig. ||, where we have used 
the input parameters given in Table |I], except that we have set = to conform to the 
limit in which these two equations are derived. The two curves differ in the large s region 
with ours becoming negative before the kinematic end-point is actually reached. 

The normalized FB asymmetry, dA{s)/ds, in the HQE-approach and the parton model are 
shown in Fig. ^ We find that this asymmetry is stable against leading order power corrections 
up to s < 0.6, but the corrections become increasingly large and eventually uncontrollable due 
to the unphysical behaviour of the HQE-based dilepton mass distribution as s approaches s"^""^ 
(see Fig. |^). Based on these investigations, we must conclude that the HQE-based approach 
has a restrictive kinematical domain for its validity. In particular, it breaks down for the high 
dilepton invariant mass region in B ^ Xgi'^i" . 
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Figure 1: Dilepton invariant mass spectrum dB{B Xse'^e^)/ds in the parton model (dashed 
curve) and with leading power corrections calculated in the HQE approach (solid curve). The 
parameters used are given in Table |I]. 
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Parameter 


Value 


mw 


80.26 (GeV) 


mz 


91.19 (GeV) 


s\v? 9w 


0.2325 


rris 


0.2 (GeV) 


rric 


i.4 [KjeM ) 




4.8 (GeV) 


mt 


175 ± 9 (GeV) 


/i 


511° (GeV) 


» (5) 
^^QCD 


0.2141°:°^^ (GeV) 


(^QED 


129 


ois{mz) 


0.117 ±0.005 


Bsi 


(10.4 ±0.4) % 
-0.20 (GeV^) 
±0.12 (GeV^) 


Ai 


A2 



Table 1: Values of the input parameters used in the numerical calculations of decay rates. 
Unless, otherwise specified, we use the central values. 



Coefficient 


Value 


cf 

ClO 


±0.3805 
-0.3110 
±4.1530 
-4.5461 



Table 2: Wilson coefficients used in the numerical calculations corresponding to the central 
values given in Table |I|. 

This behaviour of the dilepton mass spectrum in B ^ Xsi'^i" is not unexpected, as similar 
behaviours have been derived near the end-point of the lepton energy spectra in the decays 
B Xiui in the HQE approach [|16] . To stress these similarities, we show the power correction 
in the dilepton mass distribution as calculated in the HQE approach compared to the parton 
model through the ratio defined as: 

dB/ds{RQE) - d-B/ds(Parton Model) 



s = 



(55) 



di3/ds (Parton Model) 

The correction factor for B Xgi^i shown in Fig. ^ is qualitatively similar to 

the corresponding factor in the lepton energy spectrum in the decay B Xcii'g, given in Fig. 
6 of |jT6|. Finally, we note that we have been able to derive the power corrected rate for the 
semileptonic decays B Xcd-Vi obtained by Manohar and Wise in [|16|, taking the appropriate 
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Figure 2: Dilepton invariant mass spectrum dB{B XgC^e )/d,s with power corrections 
calculated in the HQE approach. The solid curve corresponds to our calculation and the 



dashed curve results from eq. (3.21) of FLS |Q with = 0. The other parameters are given 
in Table |. 
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Figure 3: FB asymmetry (normalized) (iA{B Xse^e~)/(ls in the parton model and with 
power corrections calculated in the HQE approach. The solid curve corresponds to the HQE 
spectrum and the dashed curve is the parton model result. The parameters used are given in 
Table |l]. 

limits of our calculations and taking into account the differences in our normalization of states 
and conventions, as shown in Appendix C. 

Finally, since the HQE-improved expression for the decay rate including the s-quark mass 
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effects is rather long, we give below the results in a numerical form: 

pHQE _ ^ ^^^^ ^ ^ (5g) 

where F'' is the parton model decay width for b s£^£^ and the coefficients depend on the 
input parameters. For the central values of the parameters given in Table ^ they have the 
values Ci = 0.501 and C2 = —7.425. This leads to a reduction in the decay width by —4.1%, 
using the values of Ai and A2 given in Table |l]. Moreover, this reduction is mostly contributed 
by the A2-dependent term. We recall that the coefficient of the Ai term above is the same as 
in the semileptonic width T{B ^ X^ii^e), but the coefficient of the A2 term above is larger 
than the corresponding coefficient (= —9/2) in the semileptonic decay width. Hence, the power 
corrections in F(i? ^ X^iui) and F(i? — Xgi'^i') are rather similar but not identical. 




Figure 4: The correction factor R^^^{s) (in percentage) as defined in eq. (^) for the dilepton 
mass spectrum dB{B — > Xsi~^i~)/ds. The parameters used are given in Table |l]. 

3 B-meson Wave function Effects in B ^ Xsi^i~ 

In this section, we present our estimates of the non-perturbative effects on the decay distri- 
butions in B —* Xsi~^i~. These effects are connected with the bound state nature of the B 
hadron and the physical threshold in the B Xgi^i^ in the final state. In order to imple- 
ment these effects on the decay distributions in B Xgi^i^ , we resort to the Gaussian Fermi 
motion model introduced in ||2^. In this model, the 5-meson consists of a 6-quark and an 
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spectator antiquark q and the four-momenta of the constituents are required to add up to the 
four-momentum of the 5-meson. In the rest frame of the S-meson the 6-quark and the spec- 
tator move back-to-back with three momenta pb = —pq = p. Energy conservation then imphes 
the equation 



which can only hold for all values of if at least one of the masses becomes momentum 
dependent. We treat the spectator quark niq as a momentum-independent parameter; the 
6-quark mass is then momentum dependent and we denote it by W{p): 



The 6-quark, whose decays determine the dynamics, is given a non-zero momentum having a 
Gaussian distribution, with the width determined by the parameter pp: 



with the normalization /q°° dpp'^(/){p) — 1. The distributions from the decay of the S-meson 
at rest are then obtained by convoluting the appropriately boosted partonic distributions with 
the Gaussian distribution. The resulting spectra and decay rates depend essentially on two 
parameters, pp, determining the non-perturbative width of the momentum distribution, and niq 
(or equivalently W{p)), which determines the height. In the Fermi motion model, the problem 
of negative probabilities encountered in the HQE approach for the high dilepton masses near 
s — > Smax is not present, which motivates us to use this model as a reasonable approximation 
of the non-perturbative effects in the entire dilepton mass range. The success of this model in 
describing the inclusive lepton energy spectra in S — > {Xc, Xujii^e and B Xg'y strengthens 
this hope. 

In the decay B Xs£'^£~, the distribution dB/ds depends on the Lorentz- invariant vari- 
able s only. So, the Lorentz boost involved in the Fermi motion model (Doppler shift) leaves 
the dilepton mass distribution invariant. However, since the 6-quark mass W{j)) is now a 
momentum-dependent quantity, this distribution is affected due to the difference {W{p) — mi,) 
(mass defect), which rescales the variable s and hence smears the dilepton distribution calcu- 
lated in the parton model. For different choices of the model parameters (pi?, ruq) corresponding 





(57) 




(58) 
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to the same effective 6-quark mass, (W) , the dilepton mass distributions should be very similar 
[pl| , which indeed is the case as we have checked numerically but do not show the resulting 
distributions here. 

The situation with the FB asymmetry (or the energy asymmetry) is, however, quite dif- 
ferent. Being an angular-dependent quantity, it is not Lorentz-invariant and is sensitive to 
both the Doppler shift and the mass defect. We give in Appendix E, the Dalitz distribution 
(Pr{B — > Xsi~^i')/dsdu in the Fermi motion model, given the partonic double distribution 
d'^r{b — > si^£^ )/dsdu in the 6-quark rest frame. These details, hopefully, will be useful in the 
analysis of data in 5 — > Xs£^£~ due to the popularity of the Fermi motion model. 

As we calculate the branching ratio for the inclusive decay B Xst^l" in terms of the 
semileptonic decay branching ratio B(B Xiui), we have to correct the normalization due 
to the variable 6-quark mass in both the decay rates. To get the decay rates in this model 
one first implements the wave function effects and then integrates the spectra. Fixing but 
varying the model parameters pp and rriq yields variable effective (momentum-dependent) b- 
quark mass (W). We recall that the decay widths for B Xs£~^i~ and B —>■ Xiui in this 
model are proportional to (W^) Hence the decay widths for both the decays individually 
are rather sensitive to (W). This dependence largely (but not exactly) cancels out in the 
branching ratio B{B —>■ Xsi^i^). Thus, varying (W) in the range (W) = 4.8 ±0.1 GeV results 
in Ar{B —>■ Xsi^£^)/r = ±10.8%. However, the change in the branching ratio itself is rather 
modest, namely AB{B Xsl^£^)/B = ±2.3%. This is rather similar to what we have obtained 
in the HQE approach. 

The theoretical uncertainties in the branching ratios for B Xgi^i" from the perturbative 
part, such as the ones from the indeterminacy in the top quark mass, the QCD scale Aqcd 
and the renormalization scale /i, have been investigated in the literature We have 

recalculated them for the indicated ranges of the parameters in Table |I]. The resulting (SD) 
branching ratios and their present uncertainties are found to be: 

B{B Xse+e-) = (8.4 ± 2.3) x 10"^ , 

B{B X,/iV") = (5.7 ± 1.2) X 10^^ , 

B{B ^ XsT+T-) = {2.Q± 0.5) X 10-\ (59) 

where in calculating the branching ratio B{B X^r+r^), we have included the r-lepton mass 
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terms in the matrix element |31]. These uncertainties, typically ±25%, are much larger than 
the wave-function-dependent uncertainties, and so the theoretical accuracy of the SD-part in 
the SM in these decays is not compromised by the non-perturbative effects. 

We show the resulting dilepton invariant mass distribution in Fig. ^ and the FB-asymmetry 
in Fig. |, where for the sake of illustration we have used the values {pp, mg) = (252, 300) (both 
in MeV), which correspond to an allowed set of parameters obtained from the analysis of the 
measured photon energy spectrum in B ^ X^'j , using the same model 0. We see that the 
dilepton mass distribution is very stable against Fermi motion effects over most part of this 
spectrum, as expected. The end-point spectrum in this model extends to the physical kinematic 
limit in B —>■ Xsi~^i~ , s™"^ = {itlb — mKY, which obtains for m{Xs) = itlk, as opposed to 
the parton model, in which s™-"^ = (777,5 ~ ^s)'^- The two thresholds can be made to coincide 
for only unrealistically high values of rrib and m^. The FB-asymmetry shows a more marked 
dependence on the model parameters, which becomes very significant in the high dilepton mass 
region. 

As the parameters of the Fermi motion model are not presently very well-determined from 
the fits of the existing data ^ , one has to vary these parameters and estimate the resulting 
dispersion on the distributions in B ^ Xgi^i". We show in Figs. 7 and 8 the dilepton mass 
distribution and the FB asymmetry, respectively, indicating also the ranges of the parameters 
{pF,Tnq). The resulting theoretical uncertainty in the distributions is found to be modest. 

4 LD contributions in B ^ Xgi^i' 

Next, we implement the effects of LD contributions in the processes B Xs£~^i~ . The issues 
involved here have been discussed recently in 0|-|32| and so we will be short in this part. The 
LD contributions due to the vector mesons J/ip and ip' and higher resonances, as well as the 

np 

(cc) continuum contribution, which we have already included in the coefficient Cg , appear in 
the {sLlpbiji'ii'-f^e) interaction term only, i.e. in the coefficient of the operator O9. This implies 
that such LD-contributions should change Cg effectively, but keep Cj and Cio unchanged. In 
principle, one has also a LD contribution in the effective coefficient Cj , this, however, has 
been discussed extensively in the context of the B X<j7 decay and estimated to be small 



43| , |4^ . The LD-contribution is negligible in Ciq. Hence, the three-coefficient fit of the data 
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Figure 5: Differential branching ratio dB/ds for B Xsi~^i~ in the SM including the next-to- 
leading order QCD corrections. The dashed curve corresponds to the parton model with the 
parameters given in Table 1 and the solid curve results from the Fermi motion model with the 
model parameters [pp, rriq) = (252, 300) MeV, yielding an effective 6-quark mass (W) = 4.85 
GeV. 



on B Xsi^i~ and B Xs'y , proposed in on the basis of the SD-contributions, can be 
carried out also including the LD-effects. In accordance with this, to incorporate the LD-effects 
in B ^ Xsi~^i~ , the function Y{s) introduced earlier is replaced by. 



Yis)^Y'{s)=Yis)+YrUS), 



(60) 



where Yres{s) is given as 



Yresis) 



E 



Vi=V(ls),...,t/'(6s) 



and 



(0) 



.(0) 



(61) 



(62) 



Here we adopt k = 2.3 for the numerical calculations [3C]. Of course, the data determines 
only the combination kC^^^ = 0.88. The relevant parameters of the charmonium resonances 
{^S, S) are given in the Particle Data Group [l^, and we have averaged the leptonic widths 
for the decay modes V i^i" for £ = e and i = fi. Note that in extrapolating the dilepton 
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Figure 6: Normalized differential FB asymmetry dA{s)/ds for B Xgi^i^ in the SM including 
the next-to-leading order QCD correction. The dashed curve corresponds to the parton model 
and the solid curve results from the Fermi motion model with the model parameters {pp, friq) = 
(252, 300) MeV, yielding an effective 6-quark mass {W) = 4.85 GeV. 

masses away from the resonance region, no extra g^-dependence is included in the 'J*{q'^)-Vi 
junction. (The g^-dependence written explicitly in eq. (|6T1) is due to the Breit-Wigner shape of 
the resonances.) This is an assumption and it may lead to an underestimate of the LD-effects 
in the low-s region. However, as the present phenomenology is not equivocal on this issue, any 
other choice at this stage would have been on a similar footing. The resulting dilepton mass 
spectrum and the FB asymmetry are shown in Fig. |^ and Fig. respectively. The two curves 
labeled SD and SD+LD include: 



Explicit 0(a<j)-improvement, calculated in the parton model |2^, |28[ . 



Non-perturbative effects related with the bound state nature of the B hadrons and the 



physical threshold in the final state in B ^ Xsi~^i , using the Fermi motion model 
with the parameters specified in the figures. 

In addition, the SD+LD case also includes the LD-effects due to the vector resonances, con- 

rr 

tributing to Cf as discussed earlier. 
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Figure 7: Differential branching ratio dB/ds for B Xsi~^i~ using tlie Fermi motion model for 
three different pairs of the model parameters {pp,mq) = (450,0) MeV (sohd curve), (310,300) 
MeV (long dashed curve), and {pF,mq) = (310,0) MeV (short dashed curve) yielding the 
effective 6-quark masses (W) = 4.76 GeV, 4.80 GeV, and 4.92 GeV, respectively. 

Finally, the parametric dependence due to the Fermi motion model is shown in Figs. and 
for the dilepton mass spectrum and the FB asymmetry, respectively, and compared with the 
case of the parton model in which case no wave-function effects are included. These figures give 
a fair estimate of the kind of uncertainties present in these distributions from non-perturbative 
effects. In particular, we draw attention to the marked dependence of the FB asymmetry to 
both the LD- (resonances) and wave function effects, which is particularly noticeable in the 
region rriu > m{ip'). The dilepton invariant mass spectrum, on the other hand, is very stable 
except at the very end of the spectrum, which is clearly different in all three cases shown. 



5 Concluding Remarks 

We have investigated the question of power corrections to the decay rates and distributions in 
the FCNC process B Xsi~^i~ in the HQE framework. Our motivation here was to check 
if indeed the entire dilepton mass spectrum in these decays is calculable in a theoretically 
controlled sense, which the existing results suggested [Q. Our calculations of this distribution 
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Figure 8: Normalized differential FB asymmetry dA{s)/ds for B Xgd-^d-' using the Fermi 
motion model for three different pairs of the model parameters {ppi ^q) = (450, 0) MeV (solid 
curve), (310,300) MeV (long dashed curve), and {pFyfnq) = (310,0) MeV (short dashed curve) 
yielding the effective 6-quark masses {W) = 4.76 GeV, 4.80 GeV, and 4.92 GeV, respectively. 



and the integrated rate do not agree with the ones obtained in FLS We have presented 
the details of our computations here, including the power corrections to the FB asymmetry 
not calculated earlier. In line with the analogous calculations for the lepton and photon energy 
spectra in radiative and semileptonic B decays, we have found that the HQE approach has a 
limited region of applicability in describing the dilepton mass spectrum in B Xsi^i~ . In the 
latter case, the use of the leading-order HQE approach in the high-s region results in unphysical 
distribution and hence can not be used for comparison with data. Excluding the high-s region, 
the power corrections to the dilepton mass spectrum and the FB asymmetry are, however, 
found to be small. The inclusive decay rate r{b si^i^) receives small power correction in 
the HQE approach, typically (—4%), which is similar to the one in the semileptonic decay width 
r(6 — s> uiui) but not identical. 

Despite progress in some sectors, the problem of incorporating non-perturbative effects in 
weak decays remains theoretically an intractable problem. In the present context, the structure 
functions Tj's entering the decay distributions in B ^ Xsi^i" are not known from first prin- 
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Figure 9: Differential branching ratio dB/ds for B — ^ Xsi~^i~ calculated in the SM using the 
next-to-leading order QCD corrections and Fermi motion effect (solid curve), and including the 
LD-contributions (dashed curve). The Fermi motion model parameters (Pp,mg) are displayed 
in the figure. Note that the height of the J/ip peak is suppressed due to the linear scale. 



ciples in QCD. We have modeled the non-perturbative effects in B ^ Xs 



using a popular 



Fermi motion model p2[, which allows to incorporate i?-hadron wave function effects and the 
correct threshold in the final states. Since this model gives a good description of the existing 
data on the lepton and photon energy spectra in B decays ^ |2^, we hope that it describes 
similar non-perturbative effects in the decay B Xsi~^i~ as well. We have estimated the 
dispersion on the theoretical predictions for the dilepton mass and the FB asymmetry resulting 
from the present uncertainty in the model parameters. This dispersion is marked in the high 
dilepton mass region in the FB asymmetry, but the dilepton mass spectrum is remarkably sta- 
ble. Hence, in the very high-s region, non-perturbative effects are important and have to be 
included in order to have a reliable comparison of the SM-based distributions and data, as and 
when they become available. 

Finally, we have incorporated the LD-effects using data in B decays and the measured 
properties of the resonances in the charmonium sector. As discussed in the literature, this 
is not sufficient to uniquely determine the dilepton distributions away from these resonances. 



27 



0.4 



0.3 



0.2 



0.1 



-0.1 



-0.2 



-0.3 



-0.4 



(310,300) SD 
310,300) SD-^LD 




10 



15 



20 



25 



Figure 10: Normalized FB asymmetry A{s) for B —* Xgi^i" calculated in the SM using the 
next-to-leading order QCD corrections and Fermi motion effect (solid curve), and including the 
LD-contributions (dashed curve). The Fermi motion model parameters {PF,mq) are displayed 
in the figure. 

In that context, we note that the vector-meson-dominance (VMD) approximation of the old 
vintage is often invoked to model the g^-dependence of the j-V junction, gviff)- This 
VMD-framework has been used to estimate the LD-effects in B ^ |Q. Theoretical 

uncertainties from these aspects in the dilepton mass distributions m B ^ Xgi^i' have been 
discussed in [^. We hope that data from HERA on the photoproduction of J/ip and ip' 
(and other resonances) can be used to implement the g^-dependence of the effective vertices to 
eliminate (or at least reduce) the present theoretical uncertainty from this source. However, 
as the HERA-data on the relevant processes 7*(g^) + p ^ J/iplip') +p are still preliminary 
and a g^-dependence has not yet quantitatively been extracted we do not attempt to 
undertake an improved treatment of the LD-effects in B ^ Xsi~^i~ here. In view of this, 
and the remaining theoretical uncertainties on the perturbative part discussed in section 3, the 
distributions shown in Figs. (|TT|) and (|T2|) have an overall uncertainty of order ±25%. 
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Appendices 



A The functions T>-^^(t' 4, s) 

In this appendix we list the expressions for the decomposition of the structure functions 
Ti{v.q, = 1, 2, 3) in terms of the functions Tl;-'\v.q, s), with j = 0, 1, 2, s, g, 6, representing 
the power corrections in 6 — > si~^i~ up to and including terms of order Mb /ml and explicitly 
keeping the s-quark mass dependence. The origin of these individual terms is explained in the 
text. The parton model contributions t/°^ are given in eqs. (A-1) - (A-3). 



1Mb 
X rrih 



[{1 - V ■ q)\cf T C,o\ 



(l + ml) (2 {v ■ qf - s{v ■ q) - s) - 2ml ^] \cf? 



+ - 
s 



V ■ q — s — ml {v ■ q) 



Re 



(A-1) 
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In the above expressions, tlie variable x is defined as x = 1 + s — 2 (f ■ g) — + ie. 



B Auxiliary functions Ei{s,u) and E2{s,u) in the Dalitz 
distribution (fB/dsdu{b si^i~) in the HQE Approach 



In this appendix we give the auxiliary functions Ei(s,u) and E2{s,u), multiplying the delta- 
function 6[u{s, rhs) — u^] and its first derivative 5'[u{s, rhg) — u^], respectively, appearing in the 
power corrected Dalitz distribution in 6 — s^'^^'~ given in eq. (|43|) in the text. 



Ei{s,u) = - {2 Ai 1 - Ami + - + 'drills + Amjs - 2m^s + 2m^s 

(1 - 2ml + ^t- 2^7i^s + As + s' 
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E2(s,u) 



^Ai (1 - ?fi^ + s) u{s,msf 



l_2ml + ml-f-e)(\cf\' + \C,,\^ 



+4 {1 — ml — ml + ml 



+8 1 - 2mf + - s - s Re{Cl'') 



I r^eff 1 2 

.-■2- -2 -2-2 I -2 I -2-2\ 

ims — s — ms +u +mu ' 



.eff N r^eS 



+Asu Re{Cl^) Cio + 8{i (1 + ml) Re{Cio) C^^ 



(B-2) 



C The decay rate r(6 si^i ) in the V — A hmit and 
comparison with the existing results 



In this appendix we compare our results for the power corrections in the decay B Xgt^C 
with the ones for the decays B —>■ Xcii'i, derived by Manohar and Wise (MW) |l^. In doing 
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this, we shall reduce the matrix element for the decay B Xsi~^i to the one encountered in 
B Xcii'g, obtained by the replacements: 

C9 = — Cio = - , 

Cr = 0, (C-2) 

This amounts to keeping only the CC V — A contribution in B ^ X^i^i^ . 

We remark that our hadronic states and the ones used by Manohar and Wise are differently 
normalized, with the two related by 



\M) = ^2Mb \M)^^ . (C-4) 
Hence, the forward scattering amplitudes are related through. 

Likewise, the matrix elements of the kinetic energy and the magnetic moment operators in 
the HQE approach are related, 

Ai = -2m,^K,, (C-6) 

3A2 = -2mb^Gb, (C-7) 

A1 + 3A2 = -2mb^ Eb = -2mb^Kb + Gb) . (C-8) 

Note further that our structure functions Tj are dimensionless, as opposed to the ones employed 
in [|l^. Thus, 



(TO-' ^ -if. (C-9) 

[T.r' = (c-io) 

frr, \MW (^3) rn TT\ 

~ "2M^- 

With these replacements, the structure functions Tj, i = 1, 2, 3 given in the text and Appendix 
A in this paper agree with those in MW up to the indicated normalization factors in the V — A 
limit. Note that the function Aq defined in eq. (3.4) of [jl6l transcribes to Aq = ni}^ x in our 
notation, with x defined in Appendix A. 
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After integrating T^,^ L^^ in the complex v ■ q plane, we have also compared the double 
differential distribution (l/rf,)dr/d?/dg^, given in eq. (5.2) of [|T^. Taking into account that our 
differential distributions are defined in terms of the variables u and s, as opposed to the variables 



s and y with y ^ 2 Ee used in [|T^ , and making the variable transformation y = u — 2v ■ q, we 
reproduce their result. 

Finally, using the correspondence (C-1) - (C-3), the differential dilepton distribution in 
b — i> si~^i~ reduces to (with = 0) 

^ = Q(l - sf{l + 2s) (2 + Ai) + (1 - 15s^ + 10s')X2^ , (C-12) 



r = rMi + ^Ai-^A2) , (C-13) 



which, on integration gives 

2 2 

where F^ is the parton model decay width. The above expression agrees with the well known 
result of Bigi et al [l^. Doing the same manipulation on the corresponding expressions by FLS 



33| , we get instead 

= Y'' (l + ^Ai + 13 A2 ) , (C-14) 



3 

where in F^"^"^ also only the V — A contributions are kept. This disagrees with our result as 
well as with the one in [jl5| . 



D Equivalence of FB Asymmetry and Energy Asymme- 
try in B ^ Xsi^i- 

In this appendix we address a peripheral issue, namely that the quantity Energy asymmetry, 
introduced in [RH], is simply related to the FB asymmetry, defined in |R3|, BH], and is not an 



independent quantity. It is easy to show that the configuration in which is scattered in 
the forward direction, measured with respect to the direction of the 5-meson momentum in 
the dilepton c.m.s., corresponds to the events in which > E^) in the i?-meson rest frame, 
where E± represents the £^-energy. To that end let us suppose that in the dilepton c.m.s., 
is scattered in the forward direction. In this frame, the four-momenta of C.^ and C.~ are given 
as 

Vt = , (D-1) 

P/I = (e,-P||,-P±) , (D-2) 
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where p\\ > is the longitudinal momentum and p± is the transverse momentum. Boosting the 
momenta with the velocity of B meson {v) takes one to the 5-meson rest frame, where and 
are given by 

1 V 



n 

p\\ 



(D-3) 
(D-4) 



implying that for forward scattered i+ in the dilepton c.m.s., one has E+ < E^ in the B meson 
rest frame. By using the definition of the FB-asymmetry in p^ , p5| , we obtain the following 
simple relation to the energy asymmetry of 



dAjs) 
ds 



where z = cos 9, 



1 d^B 

ds dz 

dAis 



dz 



d^B 

-1 ds dz 



dz , 



ds 



ds = BxA. 



(D-5) 



(D-6) 



where A = [N{E^ > E+) - N{E+ > E^)]/[N{E^ > E+) + A^(^+ > E^)] is the energy asym- 
metry defined in . Hence A of I^G] is identical to the normalized FB asymmetry A calculated 
in this paper. That the two quantities are related can also be seen by writing the Mandelstam 
variable defined previously, in the dilepton c.m.s. and the i?-meson rest frame: 



u 



—uis) COS 6* 



2{E+ - E_ 



(D-7) 



E 



Dalitz distribution (fT{B Xsi^i )/dsdu and FB asym- 
metry in the Fermi motion model 



We start with the differential decay rate (i^F^/ds dw dp, describing the decay b si~^i~ of a 
moving 6-quark having a mass W{p) and three momentum \p\ = p with a distribution 
which will be taken Gaussian 



dF 



B 



du — -TT-^ p (t>{p) 



1 



d^F, 
ds du' 



(E-1) 



dsdudp Ju'^^^ Mb ' ^ ^u'^ + AW{py s 

Here, d^F^/dsd-u' is the double differential decay rate of a b-quark at rest and can be written 
in the case of b —>■ si^£~ as. 



ds du' 



WtsVM 



3 a' 



192 7r3 W{p) 16 7r2 



Fi{s,p) + F2(s,p) u + F3{s,p) u 
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(E-2) 




Figure 13: Phase space boundaries for the u' and p integrations with fixed values of s and 
u drawn for s = 15 GeV^ and u = 8.9 GeV^. The integration region (sohd curve) is given 
by the intersection of (short dashed) and ±u{s,p) (long dashed curve). The fermi motion 
parameters used are {pF,mq) = (450,0). 

and the three functions have the following expressions, 
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which can be read off directly from eq. 
eff 



in the limit A,- = 0; i = 1,2. Note that the 



Wilson coefficient Cg also has an implicit W{p) dependence, as can be seen in the text. The 
integration limit for u' is determined through the equations 



Wmax = Min 



u';„ = Max 



u'+,u{s,p) 
u'_, -u{s,p) 



(E-6) 
(E-7) 



where 



Ew , _P_ 

Mr " M, 



4 s Mb^ + m2 
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(E-9) 



and 



= W s - {W{p) + nis) s - {W{p) - 



(E-10) 



A typical situation in the phase space is displayed in Fig. |I3[ Integration over p gives the 
double differential decay rate (Dalitz distribution) including the Fermi motion. The result is, 
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Note that the upper limit in p integration, pmax is determined such that p satisfies, 

^max(Pinax; "^j U) '^min(Pmax) ■S, It) . 

Lastly, the normalized differential FB asymmetry including the Fermi motion becomes. 
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where 



and 



s - (Mb + MxY] \s - {Mb - Mx f 



Mx = Max [mx, ttLs + m, 



■qj ' 



(E-13) 



(E-14) 



(E-15) 



with mq the spectator quark mass and mx the kaon mass. Since the calculations are being 
done for an inclusive decay B —>■ Xsi^i~ , we should have put this threshold higher, say starting 
from tuk + ^vr, but as this effects the very end of a steeply falling dilepton mass spectrum, we 
have kept the threshold in 5 — > Xsi~^i~ at m{Xs) = mx- 
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